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1. INTRODUCTION 

The quantum information community is currently engaged in a major effort to 
quantify the entanglement content of states of multipartite quantum systems [1-9] . A 
multipartite system is made up of a number of parts, which may be identified with in- 
dividual particles or with groups of particles. Entanglement (actually, Verschrdnkung 
is the name given by Schrodinger [10] to the nonlocal correlations responsible for vio- 
lations of the Bell inequalities [11]. This property has emerged as a physical resource 
[12,13], analogous to energy as a resource for useful work, which is to be drawn upon 
in schemes for quantum communication and quantum computation. 

Since entanglement represents uniquely quantal correlations, it becomes of great 
interest to elucidate the entanglement properties of the wave functions commonly 
used to describe strongly correlated quantum many-body systems in condensed- 
matter physics, hadronic physics, and quantum chemistry. The information gained 
in such a program should improve our understanding of quantum phase transitions 
occurring in these systems as well as their behavior in regions away from critical 
points. Here we shall take a first step in this direction by quantifying the entan- 
glement of correlated variational wave functions that have been developed to treat 
model systems of interacting Pauli spins localized on the sites of a regular lattice, 
i.e., spin-lattice models [14,15]. 

A pure state of a multipartite quantum system is entangled if and only if its 
state vector is non-separable, meaning that it is not the direct product of state 
vectors of the parts. In many-body language, reading "parts" as "particles," the wave 
function cannot be written as a product of single-particle wave functions of some 
basis. A mixed state, which may generally be represented by a density operator, 



is nonseparable and therefore entangled if and only if it is not decomposable as a 
mixture of product states. 

Exchange correlations in Fermi and Bose ground states do not contribute to 
entanglement as a useful resource [16-19]; accordingly, a state described by a single 
Slater determinant or permanent is to be considered separable. Thus, in examining 
the entanglement of a many-body wave function, one is in essence addressing its non- 
mean-field properties, which reflect fluctuation effects due to the presence of strong 
interactions. Any subset of the particles in a system of interacting particles in a pure 
state is necessarily in a mixed state. 

Bipartite (two-party) entanglement of pure and mixed states has received thor- 
ough study, especially for the case that the two subsystems are two-level systems 
or Pauli spins. (This is of course the case of most immediate concern for quantum 
computation, where the two-level computing elements are called qubits.) While the 
quantification of bipartite entanglement is well under control, analysis of multipar- 
tite entanglement quickly becomes a formidable problem as the number of parties 
increases beyond three. For an iV-partite quantum system, N > 2, entanglement 
is not characterized by a single quantity, but rather by a non-unique set of quanti- 
ties that grows polynomially with increasing N. Understandably, there is as yet no 
consensus on the best choice of such quantities. 

In the present work we will consider only (i) bipartite entanglement of a single 
spin with the rest of the spins in the lattice and (ii) bipartite entanglement of two 
spins in the lattice environment. Accordingly, our treatment of spin-lattice models 
will involve the following standard measures of bipartite entanglement [20-22]: von 
Neumann entropy, entanglement of formation, concurrence, and localizable entangle- 
ment. Considering the transverse Ising model laid out on regular lattices (square, 
cube, tessaract) in two, three, and four dimensions, information on these quantities 
will be gathered from available results on the one- and two-site (or one- and two- 
spin) density matrices corresponding to Hartree-Jastrow ground-state wave functions 
[23-25,15]. Where possible, comparison will be made with results of earlier work on 
exactly soluble models [26,27,22] or stochastic simulation methods [28]. We also make 
an interesting simple connection of the Hartree-Jastrow functions with the nilpotent 
polynomial representation of entanglement [9] , which permits us to expose important 
qualitative features of these trial ground states. 

Section II provides the necessary formal and conceptual background on entan- 
glement measures and their possible role in identifying and characterizing quan- 
tum phase transitions. In Section III we introduce the transverse Ising model and 
sketch its analysis and treatment within the framework of correlated-basis theory 
and hypernetted-chain techniques. The numerical results and attendant discussion 
are presented in Section IV. In Section V we look ahead to more ambitious explo- 
rations of the entanglement properties of the many-body wave functions employed 
in correlated-basis and coupled-cluster approaches to strongly correlated quantum 
systems. 



2. ENTANGLEMENT AND QUANTUM PHASE TRANSITIONS 



The entanglement of the two parts (1,2) of a bipartite system in a pure state 
p = may be defined as the von Neumann entropy 

S = -tr(pi log 2 pi) = -tr(pj log 2 pj) (1) 

of either subsystem, where pi = tr j(p) and pj = trj(p). When either subsystem is a 
spin-1/2 system, S ranges from (not entangled) to 1 (maximally entangled). 

The entanglement between two parts of a system in a mixed state p (e.g., two 
spins within a multispin system) is not uniquely defined. One natural definition is the 
entanglement of formation, which is the minimum, over all pure-state decompositions 
of p, of the expected entanglement required to construct p from such a decomposition, 
using S as a measure of the pure-state entanglement. In general, this quantity is 
awkward to calculate. However, for the case of two spins 1/2 (2 qubits), it can 
easily be found [21] from their density matrix as a simple monotonic function of the 
concurrence C, 

Er(p) = h (i±4HE) (2) 

with 

h(x) = —x log 2 x — (1 — x) log 2 (l — x) . (3) 

Given the two-spin density matrix p^ obtained by tracing out all spins other 
than % and j, the concurrence is calculated as 

C(pij) = max[0, p x - n 2 — A*3 — A*4] , (4) 



where the pi are the eigenvalues (in decreasing order, each real and nonnegative) of 
the Hermitian matrix 

1/2 

(5) 



R 
and 



P 1/2 ~PP 1/2 



p=(a y ®a y )p*(a y ® a y ) (6) 

is the spin-flipped density matrix (py being abbreviated as p). The concurrence C 
ranges from zero for a separable state to unity for a maximally entangled state. For 
a pure state of qubits, \ip) = a\00) + b\01) + c| 10) + c?j 1 1) , one obtains C = \ad — bc\, 
which is clearly a measure of the departure from a product state. 

Other measures of bipartite entanglement of mixed states have also been pro- 
posed. The entanglement of assistance is the maximal two-party (e.g., two-spin) 
entanglement that can be achieved by performing any kind of measurement on the 
other parts of a multipartite (e.g., multispin) system. In a sense, this measure lies 
at the opposite extreme from the entanglement of formation, and again it is hard to 
calculate. 

Verstraete et al. and Popp et al. [22] have proposed a similar measure that is more 
accessible. The localizable entanglement is the maximal amount of entanglement 
between two parties that can, on average, be created - localized - by performing 



only local measurements on the other parts of the system. Unlike the other bipartite 
entanglement measures mentioned, it is not (in general) determined from a knowledge 
of two-particle correlation functions alone. On the other hand, it captures more 
complicated features of the state of a multipartite system and leads to a natural 
definition of entanglement length. Moreover, although difficult to calculate in general, 
one can obtain bounds (usually tight ones) on its magnitude from the connected two- 
party correlation function. An upper bound is given by the entanglement of assistance 
as measured by its concurrence Ca, while a lower bound is provided by the maximal 
connected (or "classical") correlation function (see Ref. [22] for details). 

Our study of the entanglement properties of correlated wave functions for spin 
lattices was motivated by the work of Osterloh et al. [26] and more especially that 
of Osborne and Nielsen (ON) [27]. These authors were the first to explore possible 
connections between quantum phase transitions and entanglement. Both investiga- 
tions focused on the anisotropic XY model on a one-dimensional (ID) lattice with N 
sites occupied by Pauli spins with nearest-neighbor ferromagnetic couplings, subject 
to a transverse magnetic field. Since this model is exactly soluble using the Jordan- 
Wigner transform, it admits an incisive analysis of the behavior of entanglement 
in the vicinity of a simple quantum phase transition from paramagnetic disorder to 
ferromagnetic order. A special case, the transverse Ising model, received the most 
attention. There was much subsequent work on quantum spin chains along similar 
lines, driven by their tractability and by the equivalence of spin-f/2 with the qubit 
of quantum information theory. 

A quantum phase transition is associated with a qualitative change of the ground 
state of a quantum many-body system as some parameter (e.g, density, pressure, 
doping, coupling constant) is varied. In contrast to ordinary phase transitions driven 
by thermal fluctuations at finite temperature, quantum phase transitions are driven 
by purely quantal fluctuations and can occur at zero temperature. At the critical 
point in parameter space where the transition takes place, long-range correlations 
develop in the ground state. Osterloh et al. and ON proposed that there must 
exist an intimate relation between quantum phase transitions and entanglement, and 
that the behavior of a suitable entanglement measure should bear a signature of the 
singular behavior of the system near the critical point. Their results generally support 
this view, although some unexpected features were encountered. For example, the 
maximum of the nearest-neighbor concurrence does not occur exactly at the critical 
point, but at a slightly lower value of A. Importantly, these studies indicate that one 
cannot establish a universal connection between bipartite entanglement and quantum 
critical points, but rather that multipartite measures are necessarily involved in a 
rigorous analysis. 

Let us consider the transverse Ising model in the form studied by ON: 

JV-l 

H =-H (K^i + °3) ■ (7) 
i=o 

In this form, the constants defining the model are lumped into a single coupling 
parameter A. ON examined the entanglement properties of both the ground state at 
zero temperature and the thermal mixed state at finite temperature T, observing that 
the ground state has a two-fold degeneracy which is generally broken. We shall be 



concerned only with their results for the ground state, whose bipartite entanglement 
content they measured in terms of (i) the von Neumann entropy S of the one-site 
reduced density matrix and (ii) the concurrence between two spins, calculated from 
the two-site reduced density matrix. In the first case the two parties are a single spin 
ij and the N — 1 spins making up the rest of the lattice system; in the second they 
are two spins ij, residing in a mixed state within the lattice system of the remaining 
N — 2 spins. 

Single-Site Entanglement. The one-site reduced density matrix for a spin i at an 
arbitrary site (all being equivalent by translational invariance) is 

1 3 

p % = tTi(p) = -^2q a v? , 

a=0 

where a = 1 and a = 1, 2, 3 are x, y, z, and 

q a = tr(afp) = (of) . (9) 

If the full symmetries of the Hamiltonian H are enforced, the number of terms reduces 
to just one (a = 1). However, the degeneracy of the ground state leads to broken 
symmetry with qs ^ 0. The two parameters required to specify the single-site density 
matrix in the ground state are the longitudinal (x component) and transverse (z 
component) magnetizations in either of the two degenerate states, say 

M x = (0+1(7,10+), M z = (0+\a z \0+). (10) 

Thus we have 

Pi = -(I + M x a x + M z a z ) , (11) 
and the entanglement as given by von Neumann entropy is 

S = -tr(pi log 2 , 
with < S < 1. The two eigenvalues of pi are easily found to be 

« 1>2 = \ (l ± y/M* + Ml) , 

leading to the result 

S = -Kl log 2 K\ - K 2 log 2 K 2 

= -1 [(1 + x) log 2 (l + x ) + (l-x) log 2 (l - x)} (14) 

for the von Neumann entropy of a single spin with respect to the rest of the lattice, 
where 

x 2 = M* + M 2 z . (15) 



(8) 



(12) 



(13) 



Two-Site Entanglement. For the two-site reduced density matrix, similar arguments 
lead to 

3 



Pij = tT ij(p) = 7 E ® °f ' 



a,/3=0 



with coefficients 
and to the expression 



p a0 = ti(af of p) = (afaf) , 



(16) 



(17) 



Pa 



I + M z (at + ^) + J2( a »> 



a=l 



(18) 



in terms of the transverse magnetization M z and the two-spin correlation functions 
(afo-j), (a? a?), and (afap. Knowing these ingredients from an exact solution or 
approximate many-body treatment, the concurrence C may then be determined from 
Eq. (4). 



3. VARIATIONAL THEORY OF TRANSVERSE ISING MODEL 

We now review the variational-CBF approach [23-25] to the ground state and 
elementary excitations of the transverse Ising model in D dimensions. Here, "varia- 
tional" means that a variational Ansatz is made for the ground-state wave function; 
"CBF" means "correlated basis functions," implying that both the ground-state and 
excited-state descriptors will contain nontrivial correlations beyond mean-field the- 
ory. Here we are only interested in the results for the magnetizations and correlation 
functions in the ground state required for the evaluation of the relevant bipartite 
entanglement measures. 

Written for arbitrary dimension D, the Hamiltonian is written with a more 
general parametrization than that employed in Refs. [26,27,22], namely 

N N 

n = -J2 Ayof ^ + A £(1 - at) . (19) 

M i 

The iV spins are situated on the lattice sites of a .D-dimensional hypercube. A generic 
vector from one site to another will be denoted by n. The spin-spin interaction is of 
the Ising type: 

(2D n = 

A(n) = < — 1 for nearest neighbors (20) 
I otherwise 

with Aij = A(ri — rj) = A(n). The strength of the external field is measured by the 
coupling parameter A (0 < A < oo). For the special case D = 1 (which we do not 
treat numerically), identification of A with A brings the Hamiltonian (19)-(20) into 
coincidence with the form (7) used by Osborne and Nielsen, apart from an overall 
constant factor and a constant shift of energy. 



CBF theory provides a comprehensive framework for ab initio microscopic de- 
scription of strongly interacting many-body systems [29] . In application to the trans- 
verse Ising model, one would like to achieve such a description for values of the 
coupling parameter A over its full range from (corresponding to the strong-coupling 
limit) to oo (weak-coupling limit). Gross properties to be determined include the lon- 
gitudinal magnetization M x = (erf) in the normalized ground state, the transverse 
magnetization M z = (erf), the spin-spin spatial distribution function g(n) = (crfcrj), 
and the corresponding structure function S(k), all in the ground state. Further, one 
would like to determine the ground-state energy Eq and the coupling parameter A c 
at the quantum critical point, where the system changes phase from paramagnetic to 
ferromagnetic (or vice versa). (In general one would also like to find the properties 
of the elementary excitations, including the dispersion law and magnon energies.) 

To separate the mean-field effects from the effects of dynamical correlations, it 
is convenient to introduce a modified ( "connected" ) distribution function 



G(n) = (1 - Ml)' 1 [g(n) - d n0 - (1 - 5 n0 )M x 2 ] 



(21) 



and extract the so-called spin-exchange strength from M z : 



n 12 = (l-M x 2 ) 



(22) 



In mean-field approximation, G(n) = and n±2 = 1. 

The following steps have been taken in the CBF analysis of the transverse Ising 
model, and corresponding numerical results are available [23-25]: 

(i) Expression of the ground-state energy as a closed functional of the longitudi- 
nal magnetization M x , the modified distribution function G(n), and the spin- 
exchange strength n±2- 

(ii) Construction of a variational ground state having the essential correlation struc- 
ture. 

(iii) Evaluation of the spatial distribution function and spin-exchange strength, and 
hence the energy functional, for a generic trial ground state. 

(iv) Optimization of the trial ground state - derivation and solution of Euler-Lagrange 
equations. 

(v) Evaluation of the desired gross properties and correlation measures, for the op- 
timal ground state. 

(vi) Construction of the excited states and associated energies, in Feynman approx- 
imation. 

Further steps have been envisioned but not carried out: 

(vii) Systematic improvement of the zero-temperature description, by inclusion of 
higher-spin correlations and backflow in ground-state trial function and excita- 
tion Ansatz, and/or perturbation theory in a basis of correlated states. 

(viii) Extension to finite T via correlated density matrix theory [30]. 

For the purpose of the present work, only the first two steps, (i) and (ii), require 
more explicit presentation. The expression for the energy functional, applicable to a 



generic proposal for the ground state, is given by 



^[G(n),M x ;X] = (l-M 2 x ) 



D+-^A(n)G(n) 



+ A 



1-(1-M X 2 )^ 12 . (23) 



It turns out that the spin-exchange strength is dependent on G and M x . In the 
paramagnetic phase, the order parameter M x vanishes identically and Eq/N becomes 
a functional only of G. 

In mean-field theory, G(n) = and ni2 = 1, so in this case Eq becomes a 
function of M x only: 



Eq 

~N 



(M x ) = (l-M'i)D + A 1-(1-M 2 ) 5 



(24) 



This function is minimized by M x = 1 — (A/2D) 2 ] 2 , implying a critical point at 

A c = 2D, beyond which M x = 0. The resultant optimal energy is given by Eq/N = D 
in the disordered (paramagnetic) phase (A > 2D) and by Eq/N = A(l — X/4D) in 
the ordered (ferromagnetic) phase (0 < A < 2D). 

Turning to the choice of variational wave function, much of the physics of the 
transverse Ising model can be captured by a correlated trial ground state of Hartree- 
Jastrow form: 

|#hj> = exp(M x U Mx + U)\0) , (25) 



with 



i<j 

N ^ N 



(26) 



(27) 



i<j 



The vacuum or reference state |0) is taken as a symmetric product of N single-spin 
states with spin components af = +1, i.e. 



|o> = | TT ■ ■ • t>jv - 



(28) 



For a translationally invariant system, the single-spin function wi(rj) is independent 
of the lattice site r i; while the two-body pseudopotentials u and um depend only 
on the relative distance n, where n = = — r^. In mean-field approximation, 
u(n) = tt M (n) = 0. 

In the disordered phase, the Um generator is not present in the exponential form 
defining the correlated trial ground state, since the order parameter M x vanishes 
identically. However, this term makes a vital contribution in the ordered phase, 
where it is responsible for the symmetry breaking. (Note that reflection in a mirror 
plane normal to the x-axis transforms Um to — Um and reveals a two-fold degeneracy 
of the ordered ground state (characterized by the magnetizations M x and — M x ). It 
should also be mentioned that the pseudopotential «m(h) is in fact a functional of 
the generator u(n). 



Evaluation of the energy functional E [G(n, M x ; A)] requires construction of the 
spatial distribution function G(n) and the spin-exchange strength n\i corresponding 
to the trial ground state, as functionals of the pseudopotential u(n) that generates 
the spatial correlations. This is done efficiently by exploiting a 1 — 1 mapping of the 
spin-lattice system onto a binary mixture of two boson species, made possible by the 
assumed form of the trial ground state. The two boson species are characterized by 
eigenvalues +1 and —1 of the spin operator a x and may be called particles and holes, 
respectively. The partial densities p + and p_ of particle and hole components are 
determined by the magnetization M x through p± = i(l±M x ), i.e., by the expectation 

values p± = (P^) of the projectors P^ = \{1 ± erf). The Hypernetted-Chain 
(HNC) analysis available for the Hartree-Jastrow ground state of the binary boson 
mixture [31] may then be applied to determine all the requisite quantities for the 
corresponding variational description of the transverse Ising system. 

The trial ground state is optimized by deriving and solving suitable Euler- 
Lagrange equations that determine the optimal distribution function G(n) the mag- 
netization M x . For constant A and M x , the optimal G(n) is determined, through 
its HNC connection to the pseudopotential w(n), by means of the Euler-Lagrange 
equation 5E /5u(n) = for -u(n), which leads to a paired-magnon equation [24]. 
Similarly, variation of the energy functional with respect to M x at constant A and 
fixed G(n) produces an Euler-Lagrange dE /dM 2 = equation for the optimal or- 
der parameter in the ordered phase (moot in the disordered phase) , which leads to a 
renormalized Hartree equation [24]. 

This approach yields exact results in strong and weak-coupling limits, and good 
results in between, but it cannot be expected to reproduce critical exponents without 
the inclusion of higher-spin correlations. 

4. ENTANGLEMENT IN HARTREE-JASTROW GROUND STATES 

Numerical data are available from published [23-25] and unpublished variational- 
CBF calculations in the transverse Ising model that suffice for meaningful evaluations 
of measures of corresponding bipartite entanglement properties for two-, three-, and 
four-dimensional (D = 2, 3,4) versions of the model. 

Single-Site Entropy. Eqs. (14) and (15) are used to quantify the entanglement of a 
single site with the rest of the lattice (the single-spin von Neumann entropy defined 
in Section II), using the variational-CBF inputs for M x and M z for many choices of 
the coupling parameter A. The results, plotted in Fig. 1, indicate a sharp peaking of 
S{\) at the critical values of A for the order-disorder transition from ferromagnetism 
to paramagnetism given by the many-body calculation (respectively, A c = 3.14, 5.12, 
and 7.1 for D = 2, 3, and 4). It is tempting to interpret this peaking, with the 
entanglement measure reaching a maximum at the transition, in terms of a direct 
association of entanglement with quantum critical phenomena. The same behavior 
was observed by Osborne and Nielsen for D = 1, however with a distinctly higher 
maximum value of S (0.68 in comparison with the value 0.22 we find at D = 2). In 
fact, the maxima are seen to decline systematically as D increases. This finding is in 
harmony with the general understanding that classicality increases with dimension. 
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Figure 1. Von Neumann entropy S between a single site and the remaining sites 
in a square, cubic, or hypercubic lattice (respectively for dimensions D = 2, 3, or 
4), evaluated with input data from variational-CBF calculations based on optimized 
Hartree-Jastrow ground-state trial functions. The locations of the respective critical 
coupling parameters A c are labeled with (cc). 



This suggests that in practical implementations of quantum information processing 
where entanglement is used as a resource, it is advantageous to utilize chains of 
processing units rather than arrays in higher dimensions. 

Two-Site Entanglement. We next study the entanglement between two spins residing 
in the lattice in terms of the concurrence defined in Section II, deriving information 
on these measures from the data available from the CBF-variational studies based 
on the optimized Hartree-Jastrow trial function. In applying Eq. (4), We may use 
the following formulas for the eigenvalues \ii of the matrix R of Eq. (5) in terms of 
spin-spin correlation functions, which are valid in both ordered and disordered phases 



1 l! 
4 1 



> + W»l 



^3,4 = i 



1/2 



(29) 
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Figure 2. Concurrence C between nearest neighbors in a square lattice, evaluated 
with input data from variational-CBF calculations based on an optimized Hartree- 
Jastrow ground-state trial function. The location of the critical coupling parameter 
A c is labeled with (cc). The dashed line serves merely to guide the eye. 



Concurrence results for nearest-neighbor spins in the two-dimensional model are plot- 
ted in Fig. 2 at three values of the coupling parameter A. Since the existing data on 
the correlation functions is quite limited, this figure is not very informative. Still, 
it is of interest to point out that the values obtained are lower than those obtained 
by Syljuasen [28] based on Monte Carlo simulation at finite temperatures, but the 
qualitative features are the same as those found in that work and by Osborne and 
Nielsen [27] and Osterloh et al. [26]. In our case, the peak may be closer to the 
critical A than is the case for the chain [27]. As noted in Ref. [28], one should ex- 
pect the nearest-neighbor concurrence to be smaller in higher dimension due to the 
monogamous character of entanglement - the more neighbors, the smaller the share 
of bipartite entanglement allotted to each pair. Also, we know quite explicitly from 
the results of ON and Osterloh et al. that although spin-spin correlations in the usual 
sense acquire a long-range character upon approach to the critical point, this is not 
the case for the concurrence. In particular, these authors find that C vanishes for 
site separations on the chain beyond next-nearest neighbor. In our case, since the 
scale of C is already smaller because of the higher dimensionalities considered, the 
concurrence is found to vanish for pair separations beyond nearest neighbor. 
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Figure 3. Two-point connected correlation function Qa;x(| n |) for nearest neighbors 
on a square or cubic lattice (respectively for dimensions D = 2 or 3), evaluated with 
input data from variational-CBF calculations based on Hartree-Jastrow ground-state 
trial functions. The locations of the respective critical coupling parameters A c are 
labeled with (cc). The dashed line serves merely to guide the eye. 



Localizable Entanglement. As indicated in Section II, localizable entanglement is 
another recently suggested measure of bipartite entanglement [22] . It has the virtue 
of providing for a more natural definition and incisive definition of entanglement 
length than the entanglement of formation, and allowing one to observe its expected 
divergence at the quantum critical point. For the case of the transverse Ising model, 
the x-connected spin-spin correlation function Q xx (n), furnishes a tight lower bound 
on the localized entanglement, at least when nearest-neighbor spins are considered 
[22]. 

Based on the theoretical analyses and numerical data of the variational-CBF 
work reported in Refs. [23-35], we have evaluated Q xx (|n|) = <7(|n|) — M% for a 
number of A values, in the cases D = 2 and 3 (see Fig. 3). (Here, |n| denotes 
the spin-spin distance in the lattice at constant transverse magnetic field.) It is 
important to add that the data used are not fully optimal. If we compare with the 
corresponding results of Verstraete et al. [22] , it is seen that the variational results for 
Q xx behave in the same way as when, in their treatment, a small perturbing magnetic 
field is imposed in the x direction to break the parity symmetry (i.e., under x — > —x) 
and achieve a more realistic description in the strong-coupling limit (small A). The 



treatment on which our input data is based already incorporates parity breaking in 
the ferromagnetic phase. 

Tanglemeter. There have been some initial efforts toward analysis of multipartite 
entanglement in quantum spin models, focusing on spin chains [32-36]. In this vein, 
but not restricted to one dimension, we shall make some observations on the entan- 
glement properties of Hartree-Jastrow ground states within the recently developed 
multipartite formalism based on nilpotent polynomials [9]. 

Given the ground state Eqs. (25)-(27) determined by the pseudopotentials U and 
Um, one can easily find the corresponding nilpotent polynomial F defined in Ref. [9]. 
For the present argument it is irrelevant whether or not U and Um are optimized. 
It is very important that U and Um are respectively of two-body and one-body 
structure. Thus, the local contributions from Um can be ignored, since they do not 
affect the entanglement properties of the state. The required construction leading 
to the so-called tanglemeter is further simplified by the fact that the different terms 
and factors in U commute with one another. The polynomial F is then found to be 

F = 1 + ^\ /''•>< °, : + ^2 i<j<k< i ai 3 kl<T i' <r i <r k <T t + ■■■ + ai2...7va 1 +(j+ . . . a+ , 

which does not contain terms linear or of order N — 1 in the a + operators. Therefore 
F is already in its canonical form F c under general local transformations, and the 
number of the coefficients is less than the minimum number of parameters sufficient 
to describe the entanglement properties of the state. The last step is to determine 
the logarithm of the polynomial F, the so-called nilpotential /. The result for / is 
of the same form as F and so already in canonical and hence tanglemeter. 

From this simple analysis, we can infer that for any non-zero choice for the 
pseudopotential U in the form (26), the ground state (25) belongs to the general 
orbit of states which contain A-partite entanglement. However, this state cannot 
be the GHZ state containing maximal A-partite entanglement. Another implication 
is that there are no subclusters of spins that are not entangled with the rest of the 
lattice. All spins share quantum correlations with all others. 

5. ENTANGLEMENT IN OTHER MANY-BODY WAVE FUNCTIONS 

Although much can still be learned from studies of entanglement in lattice spin 
systems, the characterization of this fascinating nonlocal quantal correlation in sys- 
tems of particles having continuous spatial degrees of freedom as well as discrete 
spin/isospin observables presents the greatest opportunities for new insights relevant 
to strongly interacting quantum systems. 

A typical CBF state has the form of a correlation factor applied to a reference 
independent-particle model state. The example of the transverse Ising model shows 
that a Jastrow factor introduces nontrivial entanglement properties. It is instructive 
to compare the simple Jastrow ground state for a normal Fermi system, namely 
rii<i f(ij)® where $ is a Slater determinant, with the BCS ground state projected 
on the subspace with the same even number N of particles, which has the form 
A (x(12)x(34) • • • x(N — 1,N)) (antisymmetrized by A). The BCS wave function 



is evidently separable, being a mean-field description of Cooper pairs all in the two- 
body state x- Thus the useful entanglement is concentrated entirely in the individual 
Cooper pairs, which are not entangled one with the other. As demonstrated by Shi 
[19], the entanglement of each member of a Cooper pair with the other is given by a 
von Neumann entropy 

S = \v k \ 2 log \v k \ 2 - (1 - M 2 ) log(l - M 2 ) , (30) 

where v k is the amplitude appearing in the usual expression 

|BCS> = IJ(«* + «*4TO t -fc;)|0> ( 31 ) 

k 

for the BCS state vector, with |i> fc | 2 = (1/2)(1 - e k /E k ) and E k = (e fc + A|) 1/2 in the 
usual notation. The pair-member entanglement is zero when the gap A k vanishes 
(normal state), and the Cooper pairs are maximally entangled when E k » e k . 

By contrast, within the Jastrow product of pair functions /, entanglement is 
spread among all N(N — l)/2 pairs ij, both directly and indirectly. Inclusion of in- 
terparticle correlations of all orders (between pairs, triples, quartets, etc.) is accom- 
plished with the Feenberg function exp J2i<j u 2(ij)/2 + Y,i<j <k u 3(ijk)/2 H j, 

which in itself can provide an exact representation of the ground state of a system of 
indistinguishable bosons [37]. A study of the multipartite entanglement properties 
of this function should prove very informative. (Such an analysis may be naturally 
extended to shadow wave functions [38], a generalization of the Jastrow- Feenberg 
form used predominantly in the study of quantum fluids and solids.) 

For similar reasons, the exponential structure of the states of coupled-cluster 
theory [39], 

\y C c)=exp[S l + S 2 + ---S N ] , (32) 

wherein 

Si = ^ Cjj a|aj , S 2 = ^2c ijk iala^a k ai , etc. (33) 

ij ijkl 

is a propitious format for the investigation of multipartite entanglement in many- 
body problems in nuclear physics and related areas, already begun by Emary [40]. 
Moreover, the CC structure suggests that a fruitful analysis can be performed in terms 
of the "tanglemeter" developed by Mandilara et al. [9] as an extensive characterization 
of entanglement based on nilpotent polynomials. 
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